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We introduce an approach for quantum computing based on the theory of dynamic invariants. 
This approach generalizes adiabatic quantum computation to a nonadiabatic regime, recovering it 
as a particular case. We show that the relaxation of adiabaticity can be achieved by processing 
information in the eigenlevels of a time dependent observable, namely, the dynamic invariant oper- 
ator. Moreover, we derive the conditions for which the computation can be implemented by time 
independent as well as by adiabatically varying Hamiltonians. We illustrate our results by providing 
the implementation of both Deutsch-Jozsa and Grover algorithms via dynamic invariants. 
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I. INTRODUCTION 

Quantum information processing can be implemented 
through different quantum computation (QC) models. 
One promising such a model is provided by adiabatic QC 
(AQC) Q. In AQC, rather than using a circuit of unitary 
quantum gates as in the standard QC model (SQC), an 
algorithm is implemented via the slow continuous evolu- 
tion of a time dependent Hamiltonian H(t). The quan- 
tum system is prepared in some simple eigenstate |fi(0)) 
of the initial Hamiltonian H{Q) and is then allowed to 
evolve adiabatically so that it remains in the correspond- 
ing instantaneous eigenstate \n{t)) of H{t) at all times. 
At the end of the process, the solution of the problem is 
encoded in the final state of the system, whence it can be 
read out by means of a convenient measurement. AQC 
schemes have been proposed based on superconducting 
flux qubits [1, 11 and experimentally implemented using 
nuclear magnetic resonance techniques [4] . Protection of 
AQC against decoherence has been investigated in sev- 
eral works I, i, 0, i, i [3 , settling AQC as a favorable 
approach for QC in real (open) quantum systems. 

Inspired by AQC, the aim of this work is to propose 
a new approach to perform QC via continuous evolution 
in Hilbert space which is based on the theory of dynamic 
invariants introduced by Lewis and Riesenfeld [ll|) fl^ - 
The theory of dynamic invariants was conceived as a 
tool to solve time dependent problems in quantum me- 
chanics. In turn, as a first application, it was used to 
discuss the nonadiabatic dynamics of a time dependent 
harmonic oscillator Since then, the dynamic invari- 
ants technique has been applied to a number of problems 
which range from quantum optics [isl 0, [l^ to atomic 
systems [3) 113] • More recently, it has also been ap- 



plied to the generalization of the concept of geometric 
phases to nonadiabatic evolutions in open quantum sys- 
tems [l^ [13 ■ t'^^ present work, we will show that 
dynamic invariants can be used to implement a nona- 
diabatic approach to perform QC. In QC by dynamic 
invariants (QCDI), the computation process will be de- 
veloped in an arbitrary eigenstate (here chosen as the 
ground state) of a time dependent quantum observable 
- the so-called dynamic invariant operator, which will 
conveniently be defined below. The procedure is some- 
what similar to AQC, but it is independent of the adi- 
abatic approximation. Indeed, this is a remarkable fea- 
ture since the standard form of the adiabatic theorem [23| 
has recently been subject of controversy, with a number 
of improvements and new derivations proposed (e.g., see 
Refs. mil,!!!,!!!,!!!!). Moreover, our method allows 
for achieving the final (target) state in a nonprobabilistic 
way. This is in contrast with AQC, where the probabil- 
ity of finding out the solution is directly related with the 
adiabatic time condition 



max 

0<t<T 



{n{t)\^\m{t)) 



[En{t) - E„,{t)Y 



«1 (m^n), (1) 
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where T is the total time of evolution, H is the Hamilto- 
nian of the system, and {|'^(i))} and {i?„(i)} denote the 
set of eigenstates and eigenvalues of H, respectively p^ . 
As a by-product of our approach, we will also show how 
to recover SQC (via the simulation of constant Hamil- 
tonians) and AQC as particular cases. This is achieved 
by setting up the conditions for which the Hamiltonian 
used to implement the quantum algorithm can be chosen 
to be time independent or adiabatically varying. Hence, 
QCDI will be proposed as a method which enables us 
to bring QC to a broader scenery, regaining some other 
QC models as particular cases. As an illustration, we 
will provide implementations of QCDI for the Deutsch- 
Jozsa as well as the Grover algorithms. In particular, 
the Deutsch-Jozsa algorithm will be studied through the 
application of different kinds of interpolations between 
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the initial and final Hamiltonians, which will reveal the 
versatility of dynamic invariants for the implementation 
of QC by continuos evolution. 

The paper is organized as follows. In Section [TTl we 
review the theory of dynamic invariants. The QCDI 
method is described in Section ITIH followed by the anal- 
ysis of interpolations schemes in Section IIVI Section fVl 
illustrates QCDI in both the Deutsch-Jozsa and Grover 
algorithms. In Section IVIl conclusions and perspectives 
are briefly discussed. 



Therefore, if we initially prepare the system in the eigen- 
state \(pj{0)) of I{t) then the system will necessarily 
evolve to \ipj{t)) for any time t. The nontransitional evo- 
lution of an eigenstate of I{t) plays the role of the adia- 
batic evolution of an eigenstate of H{t) [l^, [2^. As we 
will show, this suitably built evolution in Hilbert space 
can be used to perform QC with no adiabatic constraint. 

III. QUANTUM COMPUTATION BY 
DYNAMIC INVARIANTS 



II. THEORY OF DYNAMIC INVARIANTS 

For a closed system, a dynamic invariant I(t) is defined 
as an Hcrmitian operator / which satisfies [ill, 113 



dl{t) 



[H{t),I{t)]=0, 



(2) 



dt h 

where H is the Hamiltonian of the system and, from now 
on, H will be set to one. Dynamic invariants are quantum 
mechanical constants of motion, implying therefore that 
their expectation values are constant, i.e., d{I{t))/dt — 0. 
The construction of such an operator allows for the direct 
integration of the time-dependent Schrodinger equation 



(3) 



with the dot symbol denoting time derivative. Let us 
consider an instantaneous orthonormal eigenbasis for I{t) 



I{t)\^,{t)) ^X,\^,{t)), 



(4) 



where we assume, for simplicity, that the I{t) has non- 
degenerate eigenlevels. Then, we expand the wave func- 
tion in the invariant operator basis {\ipi{t))}, yield- 
ing 



|V(i))=5]c,(t)|^,W). 



(5) 



By inserting Eq. ([5]) in Eq. ^ and projecting the result 
onto {ipj{t)\ we obtain 



(6) 



On the other hand, taking the derivative of Eq. (U) and 
projecting it onto {ipj{t)\ we get 



XiSij 



{ipj\(Pi)) = 0. 



(A, - Xj) {i{ipj\H\ip, 

The equation above implies that 

i = j : Xi — ^ Xi = constant, 

i^j- ifjl'Pi) = -i{Vj\H\ipt). 

Eq. ([7]) is a direct consequence of I{t) being a constant 
of motion. Concerning Eq. ([5]), it allows for the integra- 
tion of Schrodinger equation. Indeed, use of Eq. ^ into 
Eq. (O yields 



(7) 
(8) 



Cj(i) = Cj(0)exp 



- / dr 





d 



Let us now introduce a mechanism to perform QCDI. 
First, before the definition of the Hamiltonian opera- 
tor 7J, we introduce a time dependent dynamic invari- 
ant /(s), where s denotes the normalized time, namely, 
s ~ t/T, with T standing for the total evolution time 
and < s < 1. The operator I{s) is constructed such 
that: (a) /(O) has a nondegenerate ground state exhibit- 
ing a simple structure, (b) /(I) has a nondegenerate 
ground state which contains the solution of the problem. 
Similarly to AQC, this can be obtained by providing an 
eigenvalue penalty for any state that violates the solution 
to be found, (c) /(s) is defined, for intermediary values of 
s (0 < s < 1), by a conveniently chosen interpolation. As 
a second step, we can determine the Hamiltonian under 
which the system will be evolved by imposing that /(s) 
is a dynamic invariant. This is done here after the def- 
inition of /(s) and can be achieved by imposing Eq. ([^ 
which, in terms of the normalized time s, becomes 



dljs) 
ds 



+ iT[H{s,T),I{s)]^0. 



(9) 



Finally, we prepare the system in the ground state 1 950(0)) 
of /(O) and let it evolve during a fixed evolution time T. 
The system will then be naturally led to the correspond- 
ing ground state |iy9o(l)) of /(I), since I{s) is built (by 
definition) as a dynamic invariant. As the solution of 
the problem is encoded in |iy9o(l)), then it can be read 
out from a suitable measurement. It is worth mention- 
ing that in AQC the probability to obtain the correct 
final state at the time s = 1 goes as the inverse of the 
gap squared. Therefore, there will be in general a small 
but nonvanishing probability for AQC to fail. On the 
other hand in the nonadiabatic approach introduced here 
(QCDI) the correct final state |iy9o(l)) is reached with ab- 
solute certainty. 

The form of Eq. ^ gives a unitary evolution for /(s). 



I.e. 



/(s) = C/(s)/(0)C/t(s), 



(10) 



where U (s) is the unitary evolution operator, provided 
that H(s,T) is Hermitian, and /(O) is the invariant op- 
erator at s = 0, which one may be written as 

/(0)-l-|^o(0))(^o(0)|. (11) 

The unitary interpolation for /(s) given in Eq. (jlOp 
has the desirable feature of preserving the spectral gaps 
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among the eigenvalues of 7(s) during all the evolu- 
tion [23] ■ This ensures the absence of level crossings in 
the spectrum of I{s). Note that the ground state of /(O) 
is \(po{Q)). Then, unless a global phase factor, the final 
state of the computation |(/3o(l)) can be obtained from 
the initial state |(^o(0)) through a unitary transformation 
[/': |(^o(l)) = C/'|(/9o(0)), with U{1) = U'. Naturally, U' 
could in general be written as a sequence of simpler uni- 
tary transformations. The corresponding invariant whose 
ground state is |(/3o(l)) is given by 



(12) 



Note that no adiabaticity constraint is imposed on the 
evolution of the invariant operator /(s), and consequently 
on the state |(/9o(s)). If we allow for an adiabatic evolu- 
tion of /(s), i.e., by using that dl{s)/ds ~ 0, we obtain 
from Eq. © that the eigenstates of I{s) become also 
eigenstates of the Hamiltonian, since [ff (s, T), /(s)] ~ 0. 
Hence, we may obtain AQC as a particular case of QCDI, 
where the necessity of unitary interpolation of the invari- 
ant, given by Eq. (|10p . turns out to be unnecessary. 



IV. INTERPOLATIONS OF /(s) 

In this section we consider different kinds of interpo- 
lations of the invariant operator I{s). First we analyze 
the situation where the Hamiltonian H{s) commute with 
the evolution operator U{s) for every times s, i.e., time- 
ordering in the expansion of U{s) is unimportant. In 
this case, we will show that QC will be implemented by 
a time independent Hamiltonian. Afterwards we analyze 
the situation where H{s) does not commute with U{s). 
In this case, the theory of dynamic invariants may be an 
important tool in the design of the Hamiltonian which 
realizes the quantum algorithm. 



A. Time-ordering is unimportant: |^J/(s), (7(s) 



Let us consider an arbitrary unitary interpolation of 
the form 



U{s) — exp [—ia{s)U] , 



(13) 



where a(s) is a real function and the operator U is both 
Hermitian {U = U'^) and unitary {UU'^ = 1). By com- 
bining these two properties of U we have = 1, which 
enables us to write Eq. p3)) as 



U{s) = cos [a(s)] 1 - i sin [a{s)] U. (14) 
Substituting Eqs. ((ni) and HH) into Eq. ^U^, wc obtain 
I{s) ^ 1 ~ cos^ ais)\MO)){MO)\ 

- sin^ a(s)|(^o(l))('^o(l)| + ^ sin [2a(s)] 



For s = and s = 1 Eq. (|T5|) must be reduced to Eqs. 
(fTTj) and ([T2|) . respectively. These requirements lead to 
the following boundary conditions on a(s): 

a{0) = mr and a{l) = {n+—)Tr, with n£'Z. (16) 

Bearing in mind the conditions found in Eq. (jl6p . the 
final state |</Jo(l)) is, unless a global phase factor, given 

by 

\M^)) = u\M0))- (17) 

Note that Eq. (|16p may also be obtained from Eq. (fTil) 
by using that U{0) = 1 and U{1) |(po(0)) = U \^o{0)). To 
find the Hamiltonian that implements the interpolation 
given by Eq. (fT^ , we propose 

H{s,T) ^ h{s)U. (18) 

By using Eqs. ^ and ^ into Eq. we get 

1 da{s) 



h{s) 



T ds 



Since a{s) € R, we conclude that h{s) E R. Then, for the 
situation in which [/^ = 1 , the general form of the Hamil- 
tonian H{s, T) that generates the dynamics governed by 
operator (fT3|) is 



(19) 



Remarkably, QCDI can be implemented, under certain 
conditions, via a time independent Hamiltonian. For this 
purpose, we should choose the function a{s) as 



TTS 



with the resulting Hamiltonian reading 



H(s) = —U. 



(20) 



(21) 



X [|^o(l))(^o(0)|-|(po(0))(^o(l)|] 



(15) 



Hence we have proved that the unitary interpolation of 
the dynamic invariant allows for QC with the time in- 
dependent Hamiltonian (j2ip as long as U is unitary and 
Hermitian. By using Eq. (gTl) and Eq. into Eq. (fTU|) . 
we obtain I{t) = exp {—iHt)I{0) exp (iHt), which im- 
plies that the unitary evolution of I{t) is implemented by 
the evolution operator exp {—iHt) of the system. Note 
that the total time T appearing in Eq. (j2ip is inversely 
related to the gap of H, i.e. a larger gap will contribute to 
a smaller run time of the algorithm. This is naturally in 
agreement with AQC. Moreover, the run time in QCDI 
will be also associated with the ability of simulate efh- 
ciently the unitary transformation U, since it can involve 
nonlocal (many-body) interactions. 
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B. Time-ordering is important: 



[h{s),U(s)] /O 



We consider now a general unitary interpolation for 
the dynamic invariant which evolves |(/?o(0)) to |(po(l))- 
This interpolation can be written as 



U{s) 



exp 



N 



(22) 



with fi{s) being real functions of time and Oi time in- 
dependent Hermitian operators. Since we are concerned 
with a general form of the functions fi{s), it might be a 
hard task to obtain the Hamiltonian that implements the 
evolution operator given in Eq. (j22p . A naive procedure 
to work out H is to try a direct solution of Schrodinger 
equation 



dUjs) 



TH{s,T)U{s). 



(23) 



The application of this procedure is indeed adequate 
when time ordering of U{s) is unimportant, as shown 



above. However, for 



7^ for any time 



H{s,T),U{s) 

s £ [0, 1], the solution of Eq. ((^ may get much harder, 
since the exponential should be broken in parts composed 
by the operators Oi. An alternative route to solve this 
problem consists in the application of the theory of dy- 
namic invariants as follows. 

Let Cn = {Oi, • • • , Om} be the set of operators com- 
posing U{s) in Eq. (j^ . We assume that Cat is a subset 
of Cm = {Oi,...,Om}, with M > iV, where the ele- 
ments of Cm define an arbitrary Lie algebra 



M 
fc=l 



(24) 



where C*:- (C,f = — C^^) is the structure constant. We 



write /(O) as 



(25) 



with Ai(0) being real coefficients. Then, substituting ([22 
and ^ in (dU]), we obtain 

M ( N M 

m = E s w + * E E /^w^^^' A, (0) + ...\0k 

fc=i [ 1=1 j=i 

M 

= Y,\u{s)0,, (26) 
fc=i 

where we have used Eq. (|24p to compute the commuta- 
tors coming from the expression U{s)OkU'^ [s). As the 
operators Ok are elements of a Lie algebra, we take the 



Hamiltonian of the system as a linear combination of such 
operators 



M 



F(s,r)-E^fc(^'^)0'' 



(27) 



fc=i 



with hk{s,T) g M. This expansion of H is rather conve- 
nient since it ensures that, after evaluating the invariant 
operator (pS)) . we may obtain the coefficients hk{s,T) 
through Eq. Moreover, note that hk{s,T) can be 
determined by the solution of a set of coupled linear 
equations instead of a set of linear differential equations. 
Particularly, taking qubits as the building blocks of QC, 
we can always expand the Hamiltonian in terms of ten- 
sor product of Pauli spin matrices (satisfying the su(2) 
algebra) in the form 

H{s, T) = E ' ® ^^2^ ® ■ • • ® , 

{k^} 

where the lower index enumerate n qubits and the upper 
index refers to the set {1, ct'^, cr^, cr^} of identity and Pauli 
spin-1/2 matrices. The coefficients h'l'{s,T) g R since 
H{s,T) is Hermitian. 



V. APPLICATIONS 

As applications of QCDI, we will discuss both the 
Deutsch-Jozsa (DJ) and Grover problems. In particular, 
we will analyze the DJ algorithm for different interpo- 
lations, including the cases where H{s),U{s) = and 

H{s),U{s) 7^ 0. These illustrations will make clear 

that the theory of dynamic invariants may provide a use- 
ful tool to determine the Hamiltonian that controls the 
system dynamics during the computation process. 



A. Example 1: the Deutsch-Jozsa problem 

Given a binary function / : {0, 1}" — > {0, 1} {n is 
the number of bits) which is promised to be either con- 
stant or balanced, the DJ problem consists in determin- 
ing which type the function is. Here we construct an 
implementation by dynamic invariants for the optimized 
version of the algorithm [28j . Let us begin with the sim- 
ple case n = 1. The input state is \(po{0)) = where 
|±) = (|0) ± |l))/\/2, with {|0), |1)} being the computa- 
tional basis for the qubit (eigenstates of the Pauli matrix 
a^). The initial dynamic invariant is chosen such that 
its ground state is \(Pq{0)), i.e., 7(0) = w|— )(— |, where 
cij is a free parameter introduced to set the gap between 
the eigenstates of /(O). Note that /(O) is introduced in 
a such a way that a penalty is provided for any state 
having a contribution of |— ). Hence |(^o(0)) is its ground 
state. The DJ problem can be solved by a single com- 
putation of the function / through the unitary transfor- 
mation U\x) = (-l)-'^(^)|x) {x G {0, 1}") [13, so that in 
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the {|a;)} (computational) basis U is represented by the 
diagonal matrix U — diag[(— I)-''*-"-', (— 1)-^^^^]. In terms 
of the Pauli matrices the operator U may be written as 
t/ = ^+l+C_a^ where C± = (l/2)[(-l)/(o) ± (-1)^(1)]. 
Our implementation requires a final dynamic invariant 
/(I) such that its ground state is 



|(^o(l))=f/|¥'o(0))=e+l+>+e-|->. 



(28) 



This is accomplished by a unitary transformation on /(O), 
i.e., /(I) = UI[Q)W. Note that this is similar to the 
nonlinear interpolation for the DJ problem proposed in 
Ref. in the context of AQC (for linear AQC realiza- 
tions of the DJ algorithm, see Refs. [1^, H^])- However, 
the nonlinear interpolation is implemented here on /(s) 
instead of being realized on H{s). The final dynamic 
invariant encodes the solution of the DJ problem in its 
ground state, which can be extracted via a measurement 
of the qubits in the basis {|+), |— )}• Indeed, for a con- 
stant function, we obtain (^_|_ = ±1 and ^_ = 0. Then 
|(po(l)) = 1+) (up to a possible global phase). On the 
other hand, for a balanced function, we have ^4. = and 
^_ = ±1. Then \ipQ{l)) — |— ) (up to a possible global 
phase) . Following the reasoning developed in the Section 
HVl [see Eqs. ([TO)) and (fO)) ]. explicit evaluation of I{s) 
yields 

I{s) = - {1 - cos [2a(s)^_] 0-=" - sin [2a(s)^_] cr^} . 

^ (29) 
Since the operator U displays the properties UU'^ = 1 
and U ~ U\ we can implement the algorithm through 
the Hamiltonian 



,y(,,r) = l^t/=l^ 
T as T as 



(e+l + e-a-). (30) 



Hence, by controlling the time variation of a(s) and the 
frequency we can optimize the run time of the algo- 
rithm with no constraint in time. 



Constant Hamiltonian 



As shown in Subsection IIVAI it is possible to choose a 
time-independent Hamiltonian by adjusting a(s) = 7rs/2 
according to Eq. (PO)) . Then, the Hamiltonian becomes 



H{s,T) = —{Ul + i_a^- 



(31) 



Note that H is independent of s, but it goes as the in- 
verse of T . So QCDI for the Deutsch problem can be 
performed in an arbitrarily short time as long as we are 
able to ensure an arbitrarily large gap. Observe that H 
is diagonal. If / is constant, H is simply the identity 
operator whereas if the function is balanced, H is given 
by the Pauli matrix cr^. Note also that the Hamiltonian 
given in Eq. (pij) can be easily simulated by SQC. 



2. Oscillating Hamiltonian 

When the time derivative of a{s) is not a constant 
function, we obtain time dependent Hamiltonians. In 
order to exemplify this situation, we suppose that a{s) 
has the following form 



«(s) = fsin(f ) 



(32) 



By using Eq. (|19p. the Hamiltonian which governs the 
dynamics reads 



H{s,T)^^{s,T) (e+H-e-'T^), 



(33) 



where 



^(s,T) = ^cos(|s). (34) 

Since lo{s,T) is a time dependent function, we can com- 
pute its mean value a; (s, T) in order to obtain the effective 
frequency associated with H{s,T). Taking an average 
over the time interval < s < 1, we get 



^{s,T)^ j\{s',T)ds'^^. 



(35) 



Then, by comparing lu{s,T) with the frequency associ- 
ated with Eq. PT|) . we observe identical effective frequen- 
cies in both examples. 



3. Polynomial time dependence of the Hamiltonian 

The examples of the Deutsch algorithm explored above 
belong to the class of interpolations in which the Hamil- 
tonian commute with the evolution operator. In this 
case, despite its conceptual novelty, the practical use- 
fulness of the theory of dynamic invariants for determin- 
ing the Hamiltonian might be diminished, since H(s,T) 
can also be obtained from the evolution operator by 
directly taking its derivative. Hence, in the situation 

H{s,T),U{s) = 0, QCDI can be taken more as a 

conceptual framework for nonadiabatic QC than as a 
practical tool to derive new Hamiltonians. However, for 

H{s,T),U{s) 7^ 0, QCDI provides a systematic ap- 
proach which is inherently useful for deriving H{s, T). To 
exemplify this point, let us consider the following form 
for the evolution operator 

U{s) = exp [-^ {vs'^'a^ + /3s'V.y + 7sV^)] , (36) 

where i', 7 € M and 171,11, r G N+. In order to obtain 
the evolved state, we will rewrite U (s) as 

U{s) = exp [-id{s) ■ "ct 0(s)/2] 

= cos [0(s) /2] 1 - w{s) ■ ~a sin [(j){s) /2] , (37) 
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with 



v{s) 



l/(s) 



\v{s)\ ^l/2s2™ + /32s2n _^ ^2g2r ' 



Applying ?7(s) onto the initial state |(/?o(0)) = |+), we 
find for the final state 



{|l/(s)| cos [|l?(s)|] — ii^s™' sin 



w (s) 



(39) 



We must be able to obtain the solution of the Deutsch 
problem from |(/3o(l))- One possible way to reach this 
goal is adjusting = and /3 = 7 = (2^/2), so that 
\ipo{s)) becomes 



\ipais)) = cos 



7r^+ /s2" + s2r 



(40) 



7r^+ /s2" + s2'~ 



S" + IS' 

— I sm , 



Particularly, at time s = 1, the state of the system ac- 
quires the form 



|^o(l))=cos(^)|+> 



1 + A . /<4 

— ^ sm — — 

V2 J \ 2 



(41) 

Then, for s = 1 we obtain (unless a global phase factor) 
the state |iy9o(l)) = |+) for f{x) balanced (^+ = 0) and 
|<Po(l)) = |— ) for f{x) constant (^+ ~ ±1). In other 
words, by measuring the states |+) or |— ) we discover 
whether the function is balanced or constant, respec- 
tively. 

By taking the coefficients i', /3, and 7 as discussed 
above, the operators U{s) and I{s) become 



U{s) = exp 



'271 ^^^^ 



(42) 



and 



Us) = - 



1 — cos TT^^ 



Vs2" + s2'- 



sin 7r^+ 



Vs2" -f s2'- 



sin 7r^_f 



(43) 



To obtain the Hamiltonian responsible by the evolution 
of the state (|40p . we take it as a linear combination of 
the Pauli matrices, i.e., 

H{s, T) = K{s, T)a^ + hy{s, T)ay + h,{s, T)a\ (44) 



whose coefficients hk{s,T) {k = x,y,z) are determined 
by the replacement of Eqs. and (I44p into Hence, 
H{s,T) is found to be 



(s) = 2 |^(s)| = 2v/t^2g2m ^_ ^2^2™ _^ ^252r_ (33) ij(s, T) 



(ns2"-i -l-rs^'-i) 



2^/2T (s2" + s2r) 



2\/2(r 



Ti+r— 1 



[ns 



(45) 



Note that the Hamiltonian given by Eq. (|l5l) arises natu- 
rally from the theory of dynamic invariants while deriva- 
tion from Eq. (j23p may be clumsy. This simple example 
illustrates the advantage of dynamical invariants when 
we are dealing with general evolution operators. In the 
particular case where n = r, the Eq. (I45|) turns out to be 



His.T) 



2V2T 



(46) 



The Hamiltonian in Eq. is such that it commutes 
with U{s). For the effective frequency, we find w„(s, T) — 
7r/(2\/2r), which is of the same order as found before. 



4- Introducing more qubits: the Deutsch- Jozsa case 

We finish the discussion of the Deutsch problem by 
showing a simple implementation of QCDI for an arbi- 
trary number n of qubits. Consider the initial state 

Wm) = -i= (|0), + |1) J (|0)„ + |1)J 

where N — 2" and the state \i) is written in decimal 
basis. The invariant operator I[s) can be defined by a 
unitary interpolation as in Eq. p3p . but now with U 
defined as 



.(-1) 



/(JV-l) 



U = diag [(-1) 

N-l 

According to Eq. pT|) the invariant at s = is 

N-l 

Therefore we obtain 



(47) 



(48) 



N-l 



exp 



ITTS ^, 



(y) 



N)01, (49) 



with ^i*^'^ = (l/2)[(-l)/(') - (-I)^(J)]. The ground state 
of /(I), which will be achieved at the end of the evolution. 
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is given by 



1 



N-l 



\ipoil)) = U |^o(0)> = -^J2 (-1)^'^^ • (50) 

^ i— 

Then, by measuring the qubits in the basis {|+), |— )}, we 
wiU get that aU the qubits are |+) for a constant function 
/, otherwise the function is balanced. Concerning the 
Haniihonian, we obtain 



N-l 



2T 



(51) 



Let us analyze how efficiently we can simulate iJ as given 
by Eq. ([5T|) . Expanding H in the tensor product of the 
Pauli basis {1, cr^, cr^, ct^}, we can see that H will be 
given by a sum of diagonal terms, which are given by all 
the combinations of tensor products of 1 and cr^ for each 
qubit. These terms may involve therefore nonlocal inter- 
actions among the qubits. Since these interactions are 
in a tensor product structure we can efficiently simulate 
them [sij. However, we must guarantee that the number 
of terms in H grows only polynomially with the number 
of qubits. This is ensured due to the fact that, since the 
function is cither constant or balanced, the Hamiltonian 
will be given either by the identity operator l^r or by 
one diagonal operator of the basis given by the tensor 
products 1 and (e.g. for n = 2, H will be given in the 
balanced case by 1 (g) tr^, cr^ ® 1, or ct^ ® (t^). Therefore, 
the number of terms in H will be fixed for any n, which 
implies that QCDI can efficiently simulate H . 



Note that U is unitary and Hermitian. Let us determine 
now the conditions for which U{s) will yield an interpo- 
lation between |<i3o(0)) and the solution state \w). Indeed 
C/(0) = 1 and therefore j7(0)|^o(0)) = |<^o(0)). For the 
final time we have 



C/(l)|^o(0)) = 



{r^a - r^/3) 



\y) 



\r\ \r\ 

(55) 

Since we want t/(l)|(/3o(0)) — exp (i(/)) (where (j) is 
an arbitrary unimportant angle), we impose TxOl = rz(3. 
From this condition, we obtain 6 — e. In terms of the 
vector T*, this result implies 

Tx = P {5 + ea) , rz = a{5 + ea) \~r\ = \5 + ea| . 

(56) 

Bearing these results in mind, we are then able to build 
an oracle which allows for the determination of the ele- 
ment under search at the final time s = 1. The dynamic 
invariant I[s) can be defined by encoding |(/3o(0)) as its 
initial ground state and \w) as its final ground state. In- 
deed, this can be unitarily achieved by adopting the in- 
terpolation given by Eqs. (fTU]) and ITTI) . with U defined 
in Eq. ([54)) . yielding 



(fs) |^o(0))(^o(0)|-sin2 \w){i 



-sin(7rs)(|ii;)((^o(0)|-|^o(0))(«;|). 



(57) 



Then, the Hamiltonian operator reads 



B. Example 2: The search problem 



H 



2T 



2T 



(58) 



A simple implementation of QCDI for the search prob- 
lem [32,] can be given as follows. We start by proposing 
an oracle in a general form given by 



Uo^9\w){w\+6{\w){ipom 
+ e\ipom(MO)\, 



(52) 



so that 9,6,e G M. and \w) is the target state. The 
initial state \(po{0)) can be decomposed as \(po{0)) = 
a \w) + f3\y), where {w \y) — and + f]"^ = 1, with 
a — {w |(/Jo(0)) and P = {y |(/3o(0)) assumed as real con- 
stants. In order to rewrite Uq in terms of jw) and \y) we 
define the matrices 

1 = \w) {w\ + \y) {y\ , dy = -i {\w) {y\ - \y) {w\) , 
(^x = \w) {y\ + \y) {w\ , fj^ = \w) {w\ - \y) {y\ . 

(53) 

Then we can write Uq = tqI + r ■ a , where tq = 
(9 + e + 2aS) /2, ~a = {ax, cry, a^) and = {rx,0,rz), 
with rx ^ f3{S + ea) and = [9-e + 2a{S + ea)] /2. 
Disregarding the term proportional to the identity, we 
introduce the interpolation operator p3p with a{s) given 
by (Uni) and U by 



U = 



(54) 



Eq. dnSD into Eq. 



H 



implies 



{13(7 X 



(59) 



In order to have both a dynamic invariant and a Hamilto- 
nian which do not depend on knowing the value of\w), we 
fix a = 1/VN, which imphes (3 = y/{N- 1)/N. Then, 
Eq. ([59]) becomes 



H 



2TVN 



IcTx 



(60) 



The equation above resembles the Hamiltonian used in 
Ref. [331 to implement the analog analogue of QC. Simu- 
lation of such a kind of Hamiltonian can be implemented 
by a quantum circuit whose number of oracle calls grows 
as 0{Vn) [p. 



VI. CONCLUSION 

In conclusion, we have proposed an approach for con- 
tinuously implementing QC in Hilbert space which gener- 
alizes AQC to a nonadiabatic regime. Our method opens 



8 



up the possibility of realizing QC with new (in certain 
cases fixed) Hamiltonians and with no adiabaticity con- 
straint. Moreover it recovers SQC and AQC as particular 
cases, running in a nonprobabilistic way. For Hamiltoni- 
ans which do not commute with the evolution operator, 
QCDI was shown to provide a systematic tool to obtain 
the dynamics of the quantum system. Since QC is still 
in its early stage, the existence of a diversity of QC mod- 
els can be a valuable way of providing new experimental 
routes as well as insights for the design of quantum al- 



gorithms. Robustness of QCDI under decoherence is an 
important further point. In this context, a possible di- 
rection can be provided by the recent generalization of 
the theory of dynamic invariants for open systems evolv- 
ing under convolutionless master equations proposed in 
Ref. [l^. We leave this analysis for future research. 
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